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Abstract 

The recent measurements of lepton nucleon scattering with polarised neu- 
tron and deuteron targets are analysed together with the previous polarised 
proton data in a mutually consistent way. The detailed x-dependence of the 
polarisation asymmetry in the valence region is shown to be in agreement with 
historical predictions based on quark models. The Bjorken sum rule is shown to 
be confirmed at the Icr level and estimates of the spin content of the nucleon /S.q 
are extracted. While the average value of Ag from the three experiments comes 
out to be 0.41 lb 0.05 (to be compared with the naive quark model theoretical 
expectation of 0.58) this experimental average value is more than one standard 
deviation from the value obtained from any individual experiment. This incon- 
sistency can be overcome by allowing arbitrary higher twist contributions but 
the resulting precision is poor, Ag = 0.38 =b 0.48. 



Introduction 



The recent measurements of the polarised nucleon structure functions gi for the 
deuteron|]l| and the neutron have re-kindled the debate over the spin content of 
the nucleon which began with the measurement of gi for the proton 0] five years 
ago. The value of Jp = / g^pi{x)dx extracted in ref|0 was consistent with a tiny 
fraction of the proton's spin being carried by the constituent quarks and this fuelled 
enormous speculation over our understanding of the nucleon in the quark model 
framework. Reviews of the various interpretations of this result and of the competing 
descriptions of the proton's spin structure can be found in ref [Q. 

In order to draw conclusions from the three experiments it is important to compare 
them consistently, in particular at the same Q$2 and with the same ancillary inputs 
(e.g. the unpolarised Fi{x, Q^2) used in constructing the polarised structure function 
gi{x,Q^2) from the measured asymmetry Ai{x)). It is the purpose of the present 
paper to do this. A central plank in our analys is will be the asymmetry Ai{x) and 
we begin with a comment on this measured quantity. 

A significant feature of the data is that the x-dependences of the polari sation 
asymmetry in the valence region, Ai{x > 0.2), confirm the quark model predictions^ |^ 
for proton, neutron (see fig.l) and deuteron systems. This suggests to us that there 
is an immediate message from these data: 

The polarisation of the valence quarks is canonical 

and this should be taken into account in any attempt to interpret the data. The A{x) 
has tended to be ignored in the literature while most of the attention, and associated 
controversy, has arisen from the value of the integrated structure function gi{x) and 
its interpretation. Much of our paper will address the implications of the new data 
for this question. 

One obvious intention of a simultaneous analysis of gi^p, g^ni and g^di is to 
compare the experimental estimate of Ip^n with Bjorken's fundamental sum rule0 
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Ip.n{Q^2) = [ $lo(^i$p(a;,Q<I>2)-(7i$n(a;,Q$2)) dx = 

J 6gv 

(1) 

where we assume nj = 3. Since the three experiments carry out measurements at 
different values of Q$2 one must be careful to combine the p, n, d results a t a common 
(5$2 to test the Bjorken sum rule. For this reason and for general requirements of 
consistency we shall take only the measurements of the asymmetry from refs[0 0, ^ 
and use the latest sets of unpolarised structure functions and parton distributions to 
construct the polarised structure functions through 

. / N ^ / N Ai{x)F2{x) , . 

giix) = Ai(x)Fi(x) = — / ' ) , (2) 
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We find that g^pi (and Ip) of ref|]^ is increased as a result of tlie new information 
on F2{x) from ref|^ at low x. 

Also, in extracting 'experimental' estimates of the integrals /p,n,d we are guided by 
theoretical estimates of the asymmetry Ai at large x to cover the unmeasured region 
{x > 0.6). Even where the asymmetry is measured for x > 0.3 the experimental 
uncertainty tends to be large and can dominate the error on the integral (particularly 
SMC d data), so we prefer instead to use the valence quark model (VQM) estimates 
of Ai in this region also. 

Comparison of the Ip^n,d with the Ellis- Jaffe|P] predictions and the extraction of 
the spin content Aq of the nucleon require knowledge of the F/D parameter and 
careful treatment of QCD corrections. We re-evaluate F/D in the light of recent /3- 
decay measurements. We include the QCD corrections to the non-singlet and singlet 
contributions to the integrals. Indeed in the non-singlet case the corrections are 
known to second order at least |]10| and significantly reduce the magnitude predicted 
for the Bjorken sum rule at low Q^2. 

We find that analysing the data in this manner is consistent with the Bjorken 
sum rule at the la level. We find Aq = 0.41 ± 0.05 but the values from each of 
p, n and d lie outside the uncertainty of this mean value. Allowing for higher-twist 
contributions of arbitrary strength to force a common value of Aq from p, n and 
d leads to Aq = 0.38± 0.48. The errors on the higher twist terms themselves are 
thus large and, not surprisingly, are consistent with the rather precise theoretical 
estimates of ref[^ used in the recent analysis of Ellis and Karliner|[l^ which yields 
Aq = 0.22 ±0.10. 



Extraction of gi^p^ n, d and Ip^n,d from data 

The EMC proton experiment is at < Q^2 >~ 11 GeV$2, the SMC deuteron exper- 
iment is at < Q^2 >~ 5 GeV$2 while the SLAC E142 experiment is at < Q^2 >~ 
2 GeV<l'2. To evaluate the structure functions at a common Q^2 value of 5 GeV$2 
we take the measured values of the asymmetries Ai^p,n,d{x) for each experiment 
and assume these values hold at Q^2 = 5 GeV$2. (There is excellent evidence for 
the Q$2 independence of A^pi{x) from ref [Q] over the range 0.5 - 50 GeV$2; within 
the relatively large errors o f ref0] there is no evidence for any Q^2 dependence of 
A^di{x). Furthermore, within the precision of the SLAC E142 experiment Ai^n{x) 
appears to also to be independent of (5$2[|13|). To construct the g^ii{x,Q^2) at 
Q$2 = 5 GeV$2 we take parton distributions to compute Fi^p,n{x,Q^2) which 
are consistent with recent DIS data, in particular the F2 data of NMC IQ at small 
X. We take the Z)$/o or D^f^ distributions of MRS the latter even providing 



an excellent description of the new data from HERA[|T^. The reliable estimate of 



the gluon distribution in these fits provides, in turn, an estimate for Rqcd to insert 



2 



in eqn(^. We have checked that these results remain true when the distributions of 
ref[0 are used instead. 

The use of up-to-date unpolarised structure functions changes the values of g^pi{x) 
based on the EMC measurements at Q^2 ~ 11 GeV$2 by a significant amount; 
it is essentially the new information from NMC|§| measurements which increases 
the values of g^pi{x) and hence Ip and Aq. Computing the gi^p,n,d{x) at the 
same Q^2 value now allows us to take combinations of pairs. In fig. 2 we com- 
pare the values of xg^di{x) computed from the A^di of SMC with the combination 
^{xg^Pi{x) + xg^ni{x)) computed from the EMC and SLAC asymmetry measure- 
ments. Point-by-point we see that the two estimates are consistent with each other. 

In order to compute the integrals Ip^n,d at Q$2 = 5 GeV we must extrapolate at 
small X and large x. The small x estimate of the integral is obtained by taking the 
smallest x data point and assuming the behaviour of xgi ~ Taking a = (as 

expected from Regge behaviour) gives the central value of this extrapolation and 

the error on this is the value obtained if a = 0.5. Given that the HERA data 
on F2{x) are larger than naive Regge expectation, one may need to reevaluate the 
gi {x 0) extrapolation: our error estimate allows for some room in this direction. It 
is crucial that future experiments go to as small x as possible in order to help settle 
this question . Fig.2 indicates that the estimate for the integral / ^0.03og^di{x)dx 
from the ^{p + n) combination 

is rather different from the direct d data. 

At large x we have some theoretical guidance for the asymmetries A^p, (x) 
from valence quark models^]. Indeed, independent of the questions about the values 
of the Jj, the localised x-dependence in the valence region provides rather dramatic 
confirmation of predictions made far in advance of data on n, d even p. We regard 
this as an important clue in interpreting the polarisation data and therefore draw 
attention to, and make a brief comment on, this aspect of the data which has tended 
to receive less attention than the integral. 

As a; ^ 1 both A^pi, A^ni were predicted to reach unity ||^ 1 17 1 [?] but their values 
around x = 0.5 are expected to be quite different ||^] |Q . The expectations from the 
VQM are consistent with the measured values at x = 0.35, 0.45. We therefore use 
the VQM estimate s of Ai^p,n (with estimated uncertainties) to compute the large 
X integral (i.e. x > 0.6). In addition we notice that the final errors in Ip^n,d te nd 
to be dominated by the last two values of the gi a.t x = 0.35, 0.45 where the cros s- 
sections are small. With some caution, we choose to take the VQM values with their 
smaller uncertainties at these two x- values. Fig.2 also shows these values and we see 
that they are completely in line with the relatively precise values obtained from the 
|(p + n) combination. 

At Q$2 = 5 we compute the integrals and as a result of the above procedures 
obtain the following values:- 
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Ip{Q<^2 = 5) = 0.135 ±0.011 
/„(g<l>2 = 5) = -0.028 ± 0.006 

/d(g<l>2 = 5) = 0.041 ±0.016 (3) 

Note the value of Ip is larger than the EMC quoted value (due to new NMC 
measurements of F2) and note the larger central value of compared to that quoted 
by SMC - due to our model estimates at large x (but 1^ = 0.041 is within the quoted 
SMC uncertainty of course). From these values we can estimate in three ways the 
value of the Bjorken sum rule (where d = {p + n)/2) at Q$2 = 5GeV^2:- 

Ip-n = 0.163 ±0.013 
/2(d-n) = 0.139 ±0.035 

/2(p-d) = 0.187 ±0.040 (4) 



Bjorken Sum Rule and Spin Content of the Nucleon 

We can write for the first moments:- 

Ip = h + h + Io 
In = -h + h + Io 

Id = h + Io (5) 

where 

I3 = -^a3(l-^-3.58(^)$2) 

iZ 71 71 

36 vr TT 

/. - iaod-^) (6) 

where the next-to-leading order QCD corrections to the non-singlet quantities have 
been evaluated in ref[lC]. Note that the QCD corrections to the singlet are smaller. Q 



In eqn(y) 03 and ag are related to the F/D values while Oq is the spin fraction carried 
by quarks, i.e. 

as = F + D = Au- Ad 

as = 3F -D = Au + Ad- 2As 

ao ~ Aq = Au + Ad + As (7) 



-•^This is due to the non- vanishing of the singlet anomalous dimension 7qq$(l), S, The 0{as) 

correction to the singlet coefficient function has recently been carried outpC|]. 
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Thus we need to know F, D precisely to extract a reliable estimate for Ag. 

We have performed a fit of the current values of the /?— decay constants for np, 
Ap, SA and En|]2T| and the best fit, with x$2 = 1.55 for one degree of freedom (F/D 
with F+D constrained to equal 1.257) is 
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F/D = 0.575 ± 0.016 (8) 



F = 0.459 ±0.008 
D = 0.798^0.008 j 

This is la larger than the value used in a previous analysis of ours ||2^ principally 
due to improved Ap and Sra data. We shall use these values in what follows; however 
there are two caveats. First there is a systematic error, not included, whch arises 



from the phase space or form factor corrections in the AS* = 1 examples [ 23 1. The 
second is potentially more serious. 

The quoted figures assume that in the hadronic axial current 

icy nuQ^^lb 

= ^A7/.75 - 92 ■ (9) 

nii + rrij 

one has g2 = 0. While this is assured in the limit where rrii = rrij (such as n ^ p) 
it is not necessarily so for AS = 1. Indeed, in quark models one expects that 
92 = 0(^^^^^) with TUij = \{mi + 



'■1.3 



Hsueh et al.[^ made a fit allowing for (72 7^ and found a significant change 
in their inferred value for gA- A best fit incorporating this value raises a tantalising 
possibihty that the {9a/ 9v) t hroughout the octet are given by the naive quark model, 
all values being renormalised by 25% (such that the net spin is 0.75 rather than 1). 
Such an eventuality would correspond to the realisation of the effective quark model 
result 

F = 1/2, = 3/4 ; F/D = 2/3 (10) 

This is discussed elsewhere [p6| . These additional theoretical uncertainties merit 
further study: for the purpose of comparing most directly with the literature we shall 
adopt the 92 = best fit, eq(H). 

Since the error on F+D is so small, the uncertainty on the Bjorken sum rule 
prediction comes only from the uncertainty: we take = 0.28 ± 0.02 at Q$2 = 
5 GeV<l>2 which gives = 0.183 0.187. We note that actually the 0(a,<l>3) 
correction estimate in ref|TO] is again, like the coeffs of as and as$2, negative. In 



fig.3(a) we compare our estimates from eqn(^ with this value and we see no serious 
discrepancy. Remember that our procedure at large x produces uncertainties which 
are smaller than the true experimental errors - using the latter would further reduce 
any possible discrepancy in fig. 3 (a). 

With the above values of F/D and as, the Ellis- Jaffe0 prediction s (i.e. As = 0) 
of the integrals Ip^n,d are 
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Ip = 0.172 ±0.009 
In = -0.018 ±0.009 

Id = 0.077 ±0.009 (11) 

and these are compared with the values from eqn(^) in fig.3(b). Despite the in- 
creased estimate of Ip from the EMC data we see that only /„ is consistent with the 
assumption of As = 0. Extracting the values of Aq from eqns (^^J^J^) gives 

Jp ^ Ag = 0.21 ±0.11 (As = -0.12 ±0.04) 
/„ ^ Ag = 0.49 ± 0.06 (As = -0.03 ± 0.02) 

^ Ag = 0.24 ±0.15 (As = -0.11 ± 0.05) (12) 

to be compared with the theoretical expectation Ag = 3F — D = 0.58; these are 
shown in fig.3(c). The mean value of the spin content from the three determinations 
is 

< Ip,n,d > < Ag > = 0.41 ± 0.05 (13) 

the value being driven largely by the relatively small errors on the neutron estimate. 

Although these determinations agree at la — 2a, it is nonetheless somewhat un- 
satisfactory that the three determinations fail to give a mutually consistent value of 
the nucleon spin content. We note that the p, d, which are at moderate ly high values 
of Q$2 agree while it is the neutron data at low Q^2 that appear to be out of line. 
This calls into question the assumption made in the determination of the gi from the 
data on Ai and suggests that there may indeed be significant Q$2 dependence of 
A^ni{x) between Q$2 = 2 and 5 GeV$2 which we have neglected. We know that 
leading twist alone is insufficient to explain the unpolarised structure functio ns for 
Q^2 below ~ 4 GeV<l'2 and therefore it seems sensible to allow for some arbitrary 
higher twist contributions to the first moments of gi. Hence we write 

Ip = h + h + h + aplQ^ 
In = -h + l8 + Io + an/Q^2 

Id = Jg ± /o ± (flp ± an)/2Q$2 (14) 

We now evaluate each integral at the relevant Q^2, i.e. Ip at Q^2 = 10.7 GeV$2, 
Id at Q<l'2 = 4.6 GeV<l>2 and /„ at Q<l'2 = 2 GeV$2. Q We still use the MRS[|n 



distributions for Fi{x, Q^2) for determining gi^p, d but, as these are not valid below 

^While we are now stressing possible Q<I>2 dependence we should be concerned that even within 
the EMC and SMC experiments there is a different Q<i>2 range for each a;-value. If this is true also 
for the SLAC experiment then our higher twist analysis could be affected. 
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5 GeV$2, for the determination of gi^n{x, Q<l>2 = 2) we use both the NMC ||| and 
SLAG [27 1 data on F2 and R. As a result, we now get 



/p(Q$2 = 10.7) 
4(g$2= 2.0) 
/,(Q$2 = 4.6) 



0.134 ±0.012 
-0.023 ±0.005 
0.041 ±0.016 



(15) 



Taking a, = 0.26 at Q$2 = 10.7, a, = 0.36 at Q$2 = 2 GeV$2 we can extract the 
values of and a„ from eqn (0) by insisting on a common value of Ag from all three 
equations, using the same values of F and D as before. 
As a result we obtain 

Ag = 0.38 ±0.48 (16) 

with 



-0.161 ±0.530 
0.030 ±0.104 



(17) 



We see that the errors on Ip^n,d are such that the higher twist contributions and 
Ag cannot be pinned down with any precision. The above values in eqn(p!7D easily 
encompass the QGD sum rule estimates of Balitsky et al[|r^ used in the analysis of 
Ellis and Karliner||12||, 



= -0.005 ±0.040 
a„ = 0.039 ±0.040 



The analysis of ref[|l^ concluded that Aq = 0.22 ± 0.10 and we can see that our 
analysis indicates the sensitivity of this result to the magnitude of the higher twist 
terms. Also we see that the higher twist contribution can make a substantial reduction 
to the magnitude of the Bjorken sum rule. We note also that 0(1/Q$2) terms occur 
naturally when the Bjorken and Drell-Hearn-Gerasimovp8| sum rules are derived in 
explicit quark models[^. Indeed their magnitudes are consistent with the general 
bounds of eqn(|T7|). 

An interesting consequence of the values of Op and a„ in eqn(|l8D is that it is 
primarily the neutron which would be expected to be most affected at low Q$2 
leading to possible dramatic effects for g^ni{x,Q^2). With the above value for a„, 
/„ would be expected to change by around 50% between Q^2 of 2 and 5 GeV$2 
and we can speculate how g^ni{x,Q^2) itself would alter to bring this dramati c 
increase in the size of /„. At present the only guide we have is a compariso n of 
[xgi^d{x,Q^2 = 5) - xg^pi{x,Q^2 = 5)] extracted from SMG and EMG. There is 
a hint from this comparison that the increase would occur at very low x, especially 
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if we keep faith with the VQM prediction for x > 0.3. However since the SLAC data 
stop at X = 0.03 it is also possible that g^rii may be much larger in magnitude for 
X < 0.03, even at Q^2 = 2 GeV$2, causing In at Q$2 = 2 GeV<l>2 to be larger than 
supposed. In that case any higher twist analysis would have to be severely modified. 

Our analysis shows the importance of continued experimentation in this area. 
The situation should become clearer when SMC (who reach the smallest x values of 
all three experiments) have accumulated their full data sample and when SLAC are 
able to continue their experiment for a polarised proton target. We have stressed the 
importance 

of comparing g^p,n,di{x,Q^2) and their integrals Ip^n,d at the same value of 
Q^2 in order to pin down the spin content of the nucleon. The exercise of including 
higher twists in the analysis of present data indicates the extreme sensitivity of the 
spin content of the nucleon to their magnitude. 

In conclusion we reiterate that the successful quark model predictions of the A{x > 
0.1) imply that valence quarks are polarised canonically and that there is no need 
to rewrite the textbooks in light of these data. Immediate questions to be answered 
include whether A^n{x — 1) > [§1I|TB[; with this exception the behaviour of A{x) in 



the valence region seems establis hed and most effort is needed in the x —>■ region. In 
addition to the questions advertised above, we urge test of whether gi^d{x — > 0) < 
as this may enable a "direct" measure of Ag[p2t p6| . Finally, if the valence quarks are 
indeed polarised canonically then it becomes important to make direct measure of 
the sea polarisarion. Ref. [pO[] has argued that semi- inclusive production of fast K^— 



in polarised leptoproduction may enable the polarisation of s and/or u to be probed. 
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Figure Captions 

[1] Polarisation asymmetries A^pi{x,Q^2) and A^ni{x,Q^2) from the EMC|§[ 
and SLAC E142||^ experiments compared with the predictions of valence quark 
modeffl. The curves[|] correspond to A^pi = A^rii = with 

R = ^ and ^ = 1 (sohd), ^ = 0.75 (dashed). See refpg for further details. 

[2] xgi(^d{x,Q^2 = 5) extracted from the SMC data on y4$cii(a;) and 5fi<l>i(p + r;.)(a;, (5$2 = 
5) from the corresponding data on A^p, ni{x) from EMC and SLACE142. Also 
included are the values at large x expected from the valence quark model. 

[3(a)] Values of Bjorken sum rule extracted from the values o f Ip^n,d at Q^2 = 
5 GeV$2. The shaded region is the theoretical prediction. 

[3(b)] Ip^n,d extracted at Q^2 = 5 GeV$2 compared to the expec tation of the EUis- 
Jaffe sum rules. As = 0. 

[3(c)] Values of As and Aq = Au + Ad + As extracted from the estimates of Ip^n,d at 
at Q^2 = 5 GeV<l>2. 
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